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Abstract. We consider a functional which models an elastic body with a cavity. We show 
that if a critical point has positive second variation then it is a strict local minimizer. We also 
provide a quantitative estimate. 

S3 

1. Introduction 

The role of roughness appearing onto the surfaces and interfaces of nano-structures has been 
proved to be of great significance in several fields such as micro-electronics, metallurgy and ma- 
terials science. For instance the roughness can strongly modify the mechanical properties of mul- 
tilayered structures as confirmed by the observation that dislocations, island and cracks can be 
generated from a rough surface (see [10]). Many efforts were devoted to the investigation on how 
to control the roughness appearing onto the surfaces and interfaces of nano-structures, leading to 
the study of the so called Drive Rearrangement Instability, i.e., the morphological surfaces insta- 
bility of interfaces between solids generated by elastic stress. This phenomena has been detected, 
for instance, in hetero-epitaxial growth of thin films with a lattice mismatch between film and 
i— i substrate and in stressed elastic solids with cavities. 

The theoretical investigation of the stability of the free surface of a planar non-hydrostatically 
stressed solid has been performed in the pioneering papers by Asaro and Tiller [4] and Grinfeld [16]. 
These authors showed that the free surface is unstable with respect to a given family of sinusoidal 
fluctuation. They also gave a first insightful description of the phenomenon, nowadays named 
as Asaro- Grinfeld- Tiller instability, for which a thin film growing on a flat substrate remains flat 
up to a critical value of the thickness, after which, the free surface becomes unstable developing 
CNJ corrugations and irregularities. This instability is explained as a consequence of the presence of 

two competing energies, usually identified with a bulk elastic energy and a surface energy. After 
these results the interest of the scientific community on the rigorous mathematical study of the 
. . morphological instabilities has rapidly grown. Starting from the paper [17] where Grinfeld follows 

the Gibbs variational approach to model the morphology of thin films, it became clear that a 
second order variational analysis could be successfully used. This approach has been used in the 
context of epitaxial growth first for a one dimensional model in [6]. Then in [5] and [12] the 
model introduced in [17], which is a more realistic two-dimensional model, corresponding to three- 
dimensional configurations with planar symmetry, is studied and the problem of finding a proper 
functional setting is successfully addressed. This settled the framework in which a precise and 
detailed analysis of qualitative properties of regular equilibrium configurations has been carried 
out by Fusco and Morini in [14] via a second order variational analysis. Indeed they prove a 
sufficient condition for local minimality in terms of the positivity of second variation and provide 
a sufficiently complete picture of the phenomena that occur in epitaxially growing thin films. 

Such detailed analysis is instead far from being complete in the framework of stressed elastic 
solids with cavities. In this paper we perform a second order variational analysis for a two- 
dimensional variational model that has been recently used to describe surface instability in mor- 
phological evolution of cavities in stressed solids (see for instance [15, 19, 21]) with the aim of 
deriving new minimality conditions for equilibria and studying their stability. The model can be 
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roughly described as follows. Consider a cavity in an elastic solid, that will be identified with a 
smooth compact set Fcl ! , starshaped with respect to the origin. The solid region is assumed 
to obey to the classical law of linear elasticity, so that the bulk energy can be written in the form 



where E(u) is the symmetric gradient of the elastic displacement u and Q is a bilinear form 
depending on the material (see Section 2 for details). The surface energy is simply assumed to be 
the length of the boundary of F. Then the energy for a regular configuration is expressed by the 
functional 



In this framework the shape of the void plays a key role in the evolution of cavities in stressed 
solid bodies, while the effects of the volume changes are negligible. Hence, one usually assumes 
that the void evolves preserving its volume. The equilibria are therefore identified with minimizers 
of F(F,u) under the volume constraint \F\ = d. Since admissible configurations need not to be 
regular, the energy of such configurations has to be defined via a relaxation procedure. This 
issue, together with the study of the regularity of minima, has been addressed (even for more 
general functional involving anisotropic surface energies) in [13] where, in order to keep track of 
the possible appearance of cracks, the relaxed functional with respect to the Hausdorff convergence 
has been studied. The relaxed functional can be expressed in the following form: 



where F has finite perimeter, IV is the "regular" part of dF and represents the cracks (see 
Section 2). 

The main result of the paper is a quantitative minimality criterion that relies on the study of 
the second variation of the functional (1.1). To be more precise we prove in Theorem 5.1 that if 
(F, u) is a smooth critical configuration and the non local quadratic form d 2 J-(F, u) associated to 
the second variation of IF at (F, u) is positively defined, then there exists a constant cq such that 



for any given admissible configuration (G, v) sufficiently close to (F, u) with respect to the Haus- 
dorff distance of sets. In particular this implies not only that (F, u) is a strict local minimizer 
of (1.1) but also provides a quantitative estimate of the deviation from minimality for configura- 
tions close to (F,u) in the spirit of the recent result obtained in [1]. The minimality criterion is 
then applied to the case of a disk subjected to radial stretching where the second variation can 
be explicitly estimated to prove the local and global minimality of the round configuration if the 
applied stress is sufficiently small. 

We point out that an important open problem is how to remove the assumption of starshaped- 
ness. Indeed, even the explicit form of the relaxed functional is unknown. 

We conclude by outlining the structure of the paper and making some comments about the 
proofs. In Section 3 we calculate the second variation of F at any regular configuration (see 
Theorem 3.3) and we exploit the volume constraint to define the associated quadratic form in 
a critical configuration. At the end of the section in Lemma 3.6 we prove a "weak" coercivity 
property of d 2 F(F, u) in a critical point, which is the first step towards the proof of Theorem 5.1. 
In Section 4, as an intermediate step, we prove that the positivity of the second variation implies 
the local minimality in the C 1,1 -topology. The main point in achieving this result is to overcome 
the lack of C ' -coercivity, which would immediately imply the result. This is done by proving the 
stability of the weak coercivity under one parameter perturbation of the critical configuration (see 
Lemma 4.7). In section 5 we exploit the regularity theory for a class of obstacle problems which 
arise as perturbations of (1.1) to show that the C^-minimality actually implies the minimality 
with respect to the Hausdorff distance thus proving the theorem. In the last section we apply the 
previous analysis to the explicit case of a disk subjected to a radial stretching. 





(1.1) 




(1.2) 



F(G,v) > F(F,u)+c \GAF\ 2 
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2. PRELIMINARIES 

In this section we fix the notation and describe precisely the required background for our anal- 
ysis. We are interested in cavities identified as closed sets F with T-{}(dF) < +00 and starshaped 
with respect to the origin. The fact that F is starshaped allows us to describe it as a subgraph of 
a function. Since F has finite perimeter, the function associated to its boundary turns out to have 
bounded pointwise total variation. This will allow us to deal with functions rather than sets. 

We denote by S 1 the unit circle in R 2 and by a : K — > S 1 the local diffeomorphism defined by 
a{9) — (cos 9, sin 9), by <r _1 its local inverse and by <r ± (9) = (sin#, — cos#) its orthogonal. We set 
by C^(M) the collection of functions in C 2 (R) which are 27r-periodic. In a similar way we shall 
define the function spaces i?i(R) etc. 

With a slight abuse of notation we set 

(2.1) BV # (R) := {g : K -> [0, R ] \ g is u.s.c, 27r-periodic and pV(g, [0, 2vr]) < 00}, 

where pV(g, [0, 2ir}) is the pointwise total variation of g in [0, 27r]and i?o is the radius of a large 
ball -Br . For a function g G BV#(R) we define the extended graph of g as T g U £ 9 , where 

(2.2) T g := {pa(9) e R 2 \ g~ (9) < p < g+(9), 9 G M} 
and 



(2.3) S g = {pa{9) eR 2 \g+{8) < p< g{6), 9eR}. 

Here g~(6) := liminfg e (7(6*) and g + (9) := limsup^g 3(6*). We shall refer to S g as the sei 0/ 
cracfcs. 

Let us consider a compact set F C -Bi^ starshaped with respect to the origin. Then, for a £ S 1 , 
we can write 

F = {ra{9) eK 2 |6»e]R, 0<r < p F (9)} , 

where pf is the radial function of F and is defined by 

p F {6) := sup {p e R I pcr(6>) G F} . 

It is clear that pf '■ R — > [0, Rq] is upper semicontinuous. Moreover we have the following result, 
see [13, Lemmata 2.2 and 2.3]. 

Lemma 2.1. Let F C Br be a closed set starshaped with respect to the origin and let pp be the 
radial function of F . Then 

&F = r PF u£ PF . 

Moreover ^(dF) < +00 if and only if pp has finite pointwise total variation. 

The previous lemma rigorously shows that we may use radial functions instead of sets. Hence, 
for g G BV # {E) we set 

F g := {pa(9) eR 2 \0<p< g(8)} and tt g := B Ra \ F g . 

We may think at F g as the void and at £l g as the elastic solid. 

We can now define properly the space of admissible pairs. Given uq <E C°°(IR 2 \ Br ) we set 

(2.4) X(u ) = {(g,v) I g € BV#(R), v G H{ oc (R 2 \ F g ;R 2 ), v = u outside B B J , 

and we shall use the notation X(0) for uo = 0. We define also the following subspaces of X(uo) 

^Lip(wo) := {(9,v) G X(u ) I g is Lipschitz}, 
(2 "' 1 X Icg (u ) {(g,v) G X(«o) I .9 G C^(R), v G C 00 ^)}. 

We are now in position to give the proper definition of convergence in X(uq). 

Definition 2.2. A sequence (g n ,v n ) C X(uq) is said to converge to (g,v) in X(uo) and we write 
(g n ,v n ) — > (g,v) if 

(1) supH^dFgJ < +00, 

(2) — >• 9F 9 in Hausdorff metric, 

(3) v n — 1 « weakly in i? 1 (a;; E 2 ) for any open set w compactly contained in K 2 \ F 9 . 
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In view of [13, Lemma 2.6], we see that X(uq) is closed under the convergence of Definition 2.2. 
The elastic energy density is defined by Q(E(u)) := |C£(m) : E(u), where C is the fourth order 
tensor 

( (2 f i + A)£n + Afca 2^12 
\ 2^12 (2M+A)e22+^u 



and -E(w) is the symmetric gradient of u 



£(u) := -(Du+(Du 

The constants /1, A are called the Lame coefficients and they are assumed to satisfy the following 
cllipticity conditions 

H > and A > — /i. 

Since Q(£) > min{/z, //+ A}|£| 2 for every symmetric 2x2 matrix £, the above conditions guarantee 
that Q is coercive. We also set the ellipticity constant 

r\ := min{/i, fj, + A}. 

For a pair (g,v) G Xn p (uo) we may write the value of the functional (1.1) as 



H9,v)= I Q{E{v))dz + U l (T g 



Since this functional is not lower semicontinuous with respect to the convergence in X(uo), in 
order to effectively address the minimization problem we consider the relaxed functional 

•7 r (fl , »«) = inf{liminf J"(ff„,« n ) | (g n ,v n ) e X hip (u ), (g n ,v n ) — > (g,v)}. 

n— > 00 

The following integral representation of T is proved in [13, Theorem 3.1], where the more general 
case of anisotropic surface energy is also considered. 

Theorem 2.3. Let (g,v) € X(uo), then 

H9,v)= f Q{E(v))dz + H 1 (r g ) + 2H 1 (Z g ). 

From now on we will always deal with the relaxed functional appearing in Theorem 2.3 and 
with abuse of notation we will denote it simply by F(g, v). The minimization problem can now 
be properly stated as 

(2.6) mk{J(j,«) I (g,v) G X(u ), \Sl„\ = d} 

for some given constant d < |-Br |. Existence of solutions of the problem (2.6) is then ensured by 
[13, Theorem 3.2]. 

Given g £ BV# (R) there is one particular elastic displacement v which is the minimizer of the 
elastic energy J n Q(E(v)) dz under the boundary condition v = u outside -Br . We call this map 
the elastic equilibrium associated to g. If (h,u) G X(uq) solves (2.6) then u has to be the elastic 
equilibrium associated with h. 

Assume now that a solution (h,u) belongs to X Icg (uo) and h > 0, then (h,u) satisfy the 
Euler-Lagrange equations 

' divC(E(u)) = in O ft 

(2.7) J C(E(u))[v}=0 onT h 

^ Q(E(u)) — kh = const, on T^, 

where kh is the curvature of IV The first two equations are standard whereas the third one is the 
first variation of the functional (1.1). This motivates the following definition. 

Definition 2.4. A pair (h,u) € A rO g(u ) is said to be critical if it solves the equations (2.7). 



A QUANTITATIVE SECOND ORDER MINIMALITY CRITERION FOR CAVITIES IN ELASTIC BODIES 5 



At the end we remark that if (h, u) is a critical pair, then from the first two equations in 
(2.7) it follows that u is the elastic equilibrium associated to h. We also point out that in the 
definition of a critical point we only need to assume h to be smooth. Indeed, if we only assume 
(h,u) £ X(uo) and h G C°°(R), then it follows from the standard elliptic regularity theory (see 
[2]) that u G C°°(ft h ). 

Finally, we recall a version of the Korn's inequality which will be used throughout the paper, 
see e.g. [18]. 

Theorem 2.5 (Korn's inequality). Let ft C R 2 be a bounded domain with smooth boundary and 
v £ W ' (f2;R ). There exists a constant C = C(f2) such that if 

Dv dz — / Dv T dz, 



then 



[ \Dv\ 2 dz < C [ \E(v)\ 2 dz. 



in Jn 

Moreover if SI is an annulus A(R,r), R> r, the constant C depends only on the ratio r/R and 
C^4 as r/R^O. 

3. Calculation of the second variation 

The goal of this section is to calculate the second variation of the functional T at any point 
(h, u) £ A rcg (ito), where u is the elastic equilibrium associated to h and < h < Rq, see formula 
(3.4). We then define a quadratic form for a critical pair (3.16) and give a definition of positiveness 
of the second variation, see Definition 3.5. 

To this aim we will introduce the following notation. Given a 27r-periodic function / : K — > K 
we will denote by / : R 2 \ {0} -> K 2 the map 



(3.1) ^-'(""'(r)) 



z 

R" 

For a parameter s £ (—e,e) let (h s ,u s ) £ X reg (uo) be a smooth perturbation of (h,u), where 
u s is the elastic equilibrium associated to h s . By smooth perturbation we mean that the function 
(s, 9) i — y h s {9) is smooth and lim s _j.o ll^s - ^llc 2 (R) = 0- Moreover we set h s = ^h s , u s = and 
h' s is the derivative with respect to 9. Notations u, h means that we evaluate the time derivatives 
at s — 0. We explicitly point out that h and u depend on h s . Finally, for a given h, we define the 
set of functions 

(3.2) A(n h ) := {w : tt h -> M 2 \ (h, w) G X(0)}. 

Roughly this means that w £ A(Slh) h w — outside Br . 

We will first write the equation for ii. In the following we will denote by r the tangent unit 
vector to Th clockwise oriented and by v the unit normal to Th pointing outward the set F^. 

Proposition 3.1. Let (h,u) £ A rog (uo) be such that u is the elastic equilibrium associated to h 
and < h < Rq. Suppose (h s ,u s ) is a smooth perturbation of(h,u). Then thefunctionii £ A(flh) 
satisfies 



(3.3) 



for all w £ A(Jlh)- 

Proof. The proof is very similar to the one in [14]. Arguing as in [7, Proposition 8.1] we obtain a 
one parameter family of C°°-diffeomorphisms 3> s ( - ) : K 2 \ {0} — > K 2 \ {0} such that <& s (z) = h s on 
dF h . ~ 



CE(u) : E(w) dz = I (h, v) CE(u) : E(w) dU 1 

div T ({h,u)CE(uj) -wdH 1 , 
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Suppose first that w £ A(flh) H C°° (f2/j). We may extend w outside flh such that w € A(Slh a ) H 
C°°(r2^ s ). Since m s is the elastic equilibrium in f2^ s we have 



CB(u s ) : E{w)dz = 0. 
Differentiate this with respect to s and evaluate at s — to obtain 

/•2-7T 

CE(u) : E(w) dz- hh [CE(u) : E(w)} (h <r(6)) d8 = 0. 
Jo 

Using the area formula and notation (3.1) we may write 

CE(u) : E(w)dz = [ (h, v) CE(u) : E{w) dU 1 , 



where we have used the fact that the normal can be written in polar coordinates as v = h °~t h ° . 
The rest will follow by integration by parts and from the fact that CE(u)\v\ = on Th as in (2.7). 
To obtain (3.3) for every w G A(Clh) one may use a standard approximation argument. □ 

Remark 3.2. Notice that the equality (3.3) clearly holds also for test functions of the form 
w(z) = w(z) + Az + b, where w G Aiflh), b € M 2 and A is an antisymmetric matrix. 

In the next theorem we derive the formula for the second variation of J 7 . 

Theorem 3.3. Suppose that (h,u) and (h s ,u s ) are as in Proposition 3.1. Let v be the outer 
normal ofTh = dFh, t be the tangent (positively oriented) ofTh and k be the curvature ofTh- 
The second variation of IF at (h,u) is 

^F(h s ,u s )\ s=Q =- [ 2Q(E(u))dz+ I \d T (h,v)\ 2 dU 1 



(3.4) 



(d v Q{E(u)) + k 2 ) (h, v) 2 dU x 
(Q(E(u))-k)d T ((h,v)(h,T)) dU 1 
(Q(E(u))-k) (M>) dn\ 

Proof. We will treat the elastic and the perimeter part separately and write 

F{h s ,u s )= [ Q(E{u s ))dz + H 1 (r hs )=F 1 {h s ,u s )+F 2 {h s ). 

Since h s is smooth, we notice that S/ ls = and denote by $ s the family of diffeomorphisms as in 
the proof of Proposition 3.1. 

1st Variation : We start by differentiating the perimeter part ^(hg). 



Since H 1 ^,) = y/h 2 + h' 2 d6 we have 



d _ „ . f 2n h s h s + h'h' 



1 ^2(h s )= I ' s , s 3 d6. 
ds s> J ^Jhl + hf 



Integrate the second term by parts and obtain 



o JhJ+W 2 Jo KVW+h! 2 (ht + hff/2 



Then we have 



(3.5) 



d C 27r ■ ( h 2 + 2h' 2 -h h"\ /' 27r 



(hs, vh„) ksdH 1 . 
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where k s — ^|^2 2 ^/ i /2)3/2 a is the curvature of Th 3 in polar coordinates. 

Let us now treat the elastic part T\{h s , u s ). Differentiate it with respect to s and get, as in the 
proof of Proposition 3.1, 

^T 1 (h s ,u s )= [ CE{u s ):E{u s )dz~[ h s h s Q{E{u s )){h s a) d0. 
ds Jn hs Jo 

The first term disappears since u s satisfies the Euler-Lagrange equations (2.7) and ii s € A{flh B )- 
Hence, we are only left with 

(3.6) ^F 1 {h s ,u s ) = - I h s h s Q{E{u s )){h s a)d6 = - [ (h s , v K ) Q{E(u s )) dU 1 . 

ds J J Thg — 

Combining (3.5) and (3.6) gives the first variation of T . 
2nd Variation : We will divide the proof in two steps. 
Step 1: As in [14], we begin by making a couple of general observations. 
Let d be the signed distance function from T^, i.e., 



d(z) 



-dist(z,I\) if z e F h , 
dist(z,T /l ) Xz#F h . 



Since the boundary Th is a graph of a C°°-function, d is C 1 function in a small tubular neigh- 
bourhood of Th t - Setting v{z) := Vei(z) and k{z) := (div^)(z), we observe that on Th, v is the 
outer normal to and k is the curvature of T^. 
First we claim that 

(3.7) d v k = -k 2 , on T h . 

Differentiating the identity \v\ = 1 with respect to v yields Dv[v] = 0. This shows that 

(3.8) Dv = D T v = kr ® r and div v = div r v, onT/,. 

Differentiating the identity Dv[v] = yields Y^j = i{d 2 k ViVj + djVid k Vj) = for k,i = 1,2. Hence 
we have 

2 2 

for i, k = 1,2. Using the previous identity we obtain 

d v k = Trace {d v {Dv)) = -Trace {{Dv) 2 ) = -k 2 on Th, 

where the last equality follows from (3.8). Hence we have (3.7). 
Next we claim that 

(3.9) (v,T) = -d T {h,v), on I\. 

Recall that < & s (z) : K 2 \ {0} — > M 2 \ {0} is a one-parameter family of C°°-diffeomorphisms such 
that &s{z) = h s on Th- Notice that we have 

(3.10) (*,!/) = (h,v), onT h . 

Differentiating D^ T D^[v] — v and calculating at s = gives D$~ T [v} — — D$ T [v}. Differen- 
tiate the identity 

v s o $ s = 8 J- 

\D$7 T [v]\ 

with respect to s, evaluate at s = and use the previous identity to obtain 

(3.11) v + Dv[<5>] = -D$ T [v\ + {D^> T [v],v)v, onT h . 

By (3.8) we have Dv = D T v T on Th- Therefore, multiplying (3.11) by r we obtain 

(z>,r) = -(D$ t [v],t) - (£>!/[$], r) 
= -(.D$ t [i/],t) - (L^ t [$],t) 
= <(-£)<$, i/) ), r) =-d T {h, u) onT h 
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and (3.9) is proven. 

Step 2: Let us start with the perimeter part and differentiate (3.5) 



-7-2 ^2 (M| =n = / hhk(ha)d8+ hh 2 d a k(ha)d8 
ds s Jo Jo 

h 2 k{ha)d9 + / hhk(ha)d6. 



JO 



For the term A we have that 

[■2n 



A= [ hhk{ha)d6= [ (h,v)kdU 1 = [ (h,u) AwrvdU 1 
Jo Jr h Jv h 

(v,T)d T (h,v)dH 1 = { \d r {h,u)\ 2 dU l , 



where we have used (3.9). For the term B, noticing that 

d a k = — d v k ; d T k and r = — , , 

Vh 2 + h' 2 Vh 2 + h' 2 Vh 2 + h' 2 

we may write 

B= [ hh 2 d a k(hcr)de= [ (h,v) 2 dvkdH 1 + [ (A, v){h, r) d T kdU x 
Jo Jr h Jr h 

{h,u) 2 k 2 dU 1 - [ kd T ((h,is)(h,T}) dH\ 
Jr h v ' 

where we have used (3.7) and integration by parts. Hence, we have 
d 2 r r ■ 



/2 ^2(ML =0 =/ \dr{h,v)\ 2 dU L - (h^yk'dH 1 

(3.i2) " Jr ; hh . ,. 



;d T ({h,v){h,T)) dU 1 + [ k { ^^-dU 1 +( k(h,v)dU 1 . 



We are left with the elastic part. Differentiate (3.6) to obtain 

j2 fl-K i-2-k 



ds 2 



P Z7T PA7T 

T 1 (h s ,u s )\ s=Q = - CE(u) : E(u)hhdO- d a Q(E(u)) h h 2 d9 
Jo Jo 



2jv 



Q{E{u)) (A 2 + hh) dO. 
Since u G A(Qh), we may rewrite the first term using (3.3) as follows 

2-7T p 

CE(u) : E{u)hhdO = \ (h,v) CE(u) : E{u) dU 1 

2Q{E{ii))dz. 

JQ h 

For the second term, noticing that 

»)) - Tp^i WW) ~ 7 j^= drQ{E{u)) 



and using integration by parts, we get 

/.2tt 



/ d T Q(E(u)) 


((A,^)(A,r)) 


f Q(E(u))d T 
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Finally we have that 
d 2 



js2 Ti{h s ,u s )\ s=Q = - / 2Q(E{u))dz- I d v Q{E(u)) {h,uf dU l 



(3.13) +/ Q(E(u))d T ((h,v)(h,T)) dH 1 - Q(E(u)) [ -=^ydH 1 

- / Q(E{y))(k,v)dH}. 
Combining (3.13) with (3.12) yields the formula (3.4). 



□ 



In the formula (3.4) we considered any smooth perturbation h s of h. However, in order to be 
admissible for our minimization problem, a perturbation h s has to satisfy the volume constraint 
\Fh s \ = \F h \, or equivalently 

(3.14) / h 2 s d6 = h 2 d9 for all s > 0. 



Remark 3.4. If (h,u) £ X rcg (uo) is a critical pair and the perturbation (h s ) satisfies the volume 
constraint (3.14), then the last two terms in (3.4) vanish. Indeed one term vanishes because the 
term Q(E(u)) — k is constant on Th by (2.7). The second one vanishes since differentiating two 
times the volume constraint (3.14) with respect to s we obtain 



(k, 



2 



+ (h, v) dH 1 = . 



Motivated by the previous observation, for any i\> g (K) satisfying 
(3.15) / hi/;d6 = 0, 



Jo 

we define the quadratic form associated to a regular critical pair (h, u) 

d 2 F(h,u)[i>] :=- [ 2Q(E(u, 4 ,))dz+ f \d T {±, v) \ 2 dU 1 

(3.16) Uh h 

- / (d,Q(£(u)) + fc 2 ) (^} W , 

where € A(flh) is the unique solution to 

(3.17) / CE(u i ,):E(w)dz = - I div T v) CE{u)) ■ w dH 1 , Vw e A{fl h ) . 
Jn h Jr h 

We define now what we mean the second variation of T being positive at a critical pair. 

Definition 3.5. Suppose that (h,u) € X Icg (uo) is a critical pair. The functional (1.1) has positive 
second variation at (h, u) if 

d 2 F{h,u)[f[ > 
for all ip £ Hl(M.) such that ip ^ and satisfies (3.15). 

We point out that if the second variation is positive at a critical point (h,u), then the formula 
(3.4) and Remark 3.4 imply that for every smooth perturbation h s of h satisfying the volume 

constraint ^^(hs, u s )\ s=0 > 0. 

At the end of the section we prove the following compactness result. 

Lemma 3.6. Suppose that a critical pair (h, u) £ X Icg (u ) is a point of positive second variation, 
and < h < Rq ■ Then there exists Co > such that 

d 2 F(h,uM>co\\(±,v)\\ 2 HHrh) , 
for every tp £ H^(M) satisfying (3.15). 
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Proof. First notice that the condition (3.15) can be written using the notation (3.1) as 
(3.18) / {ijj, v) dn 1 = 0. 



Using the Sobolev-Poincare inequality IKV 7 : ^)llL 2 (r h ) < C\\d T (il>, f)||L 2 (r h ) an d (3.18) we easily 
see that it suffices to show that 

c := inf j<9 2 J"(/i,w)[V>] | tp <= satisfying (3.15), J \d T {±,v)\ 2 dU 1 = l| > 0. 

Choose a sequence (VVi) such that ?/>„ are smooth, satisfy (3.15), J r \d T (ip n ,i/)\ 2 dT-L l = 1 and 

By restricting to a subsequences, we may assume that (ip n , v) — / weakly in H 1 ^^). By denning 

(a, i/; n(v) 

we see that / = for some ^ € i?^(M). Moreover since J r / dH 1 = 0, the function -0 

satisfies (3.15). 

Next we prove that T(h, u) has the following lower semicontinuity property 
(3.19) lim d 2 F{h,u)[4> n ] > d 2 F(h, 



Indeed, since (ip n , v) — (tp, v) weakly in H x {£h) then (ip n , v) — > (0, z/) strongly in L 2 (Th). There- 
fore we only need to check the convergence of the first term in (3.16). 

First of all, the smoothness of ip n implies that u^ n is smooth. Consider a domain = B2R \Ff i 
and a map w n (z) = u^ n (z) + A n z + b n , where A n is an antisymmetric matrix and b n € K 2 . Notice 
that w n £ iJ 1 (f2/ l ). By choosing A n such that Dw n dz — f^Dw^dz we have by Korn's 
inequality (2.5) that H-D^nH^^) — C \ \E(w n )\\ L2 (Q h y Moreover, since v^, n = outside Br 01 we 
have ||£'(w n )|| i 2(f2 h ) = \\E( u il: n )\\L 2 (n h )- By the Remark 3.2 we may use w n as a test function in 
(3.17) and using Holder's inequality and the trace theorem we get 



2Q(E{v^, n )) dz = - div r ((^, v) <CE(uj) ■ w„ dH 1 

(3_ 20) < \\(tn,v)CE(u)\\ Hl{rh) \\w n \\ L2{rh} 

<C\\($n,v) CE( U )\\ H1(n jDw n \\ LHClh) 
<C\\{^,v) £E{u)\\ Hl{nh) \\E{u^)\\ L ^ h y 

Therefore 

\\Dw n \\ L , {hh) <C\\E{u^)\\ L , [nh) <C. 
However, since = outside Br we get 

\B 2Ro \ B Ro \ \A n \ 2 = f \E(w n )\ 2 dz < C. 

J B 2 R \B Rg 

This implies that the matrices A n are bounded and therefore H-Du^H^n^) < C. 

By the compactness of the trace operator we now have that u^ n — > in L 2 (Th) up to a 
subsequence. Use as a test function in (3.17) to obtain 



lim / 2Q(E(utpJ) dz = - lim / div T ((ip n , v) CE(u)) ■ iu n dU 1 

div r ({<ip,v) CE(u)) ■ u^dU 1 
2Q(E{ Ui ,))dz. 



This proves (3.19). 
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The claim now follows since if ip ^ 0, the lower semicontinuity (3.19) implies 
c = lim &F{h,u)\i) n ] > d 2 F(h,u)hp} > 0. 

n— >oo 

On the other hand if = then the constraint J T \d T {ip ni v) | 2 dH 1 — 1 yields 

c = lim d J-(h^ u)[ip n ] = 1. 

n— >oo 

□ 



In this section we perform a second order analysis of the functional (1.1) with respect to C ' - 
topology in a spirit of [11]. The main result is Proposition 4.1 where it is shown that a critical point 
(h,u) € X reg (uo) with positive second variation is a strict local minimizer in the C 1, """-topology 
and the functional satisfies a growth estimate. We point out that, according to Lemma 3.6, the 
second variation at (h,u) is coercive with respect to a norm which is weaker than the C ' -norm. 
Therefore the local minimality does not follow directly from Lemma 3.6. The idea is to prove a 
coercivity bound in a whole C 1, ^neighborhood of the critical point, which is carried out in Lemma 
4.7. The main difficulty is to control the bulk energy, that will be done by using regularity theory 
for linear elliptic systems. We prove the main result at the beginning without worrying about the 
technicalities. All the technical lemmata are proven later in the section. 

Proposition 4.1. Suppose that the critical pair (h,u) € X Tes (uo) is a point of positive second 
variation such that < h < Rq. There exists 5 > such that for any admissible pair (g, v) € X(uq) 
with g e C^ l (R), | \g\ U=([o,27r)) = IMIl 2 ([o,2tt)) and \\h - g\\ c i,i (R) < 6 we have 

H9,v) > •? 7 (M)+ci||/i-3||i2 ([0)2 , r)) . 

Proof. Assume first that g <E C^?(R) and \\h — g\\c 2 (R) < 8. By scaling we may assume that 

i 

\\h\\m[ 0< 2-«)) = (jftfdBy = 1. We define 

h + t(g - h) 
9t := \\h + t(g-h)\\ L 2 
so that gt satisfies the volume constraint, and set 

f(t) :=H9t,Vt), 

where v t are the elastic equilibria associated to g t . We calculate 4^T{gtiVt) for every t € [0, 1) 
by applying the formula (3.4) to (gt) s — 9t+s of gt and get 

cP_ 

dt 2 ./,, 

(d Vt Q{E{v t )) + k 2 ) {g_t,v t ) 2 dU 1 

t 

(Q(E(v t ))-k t )d Tt ((gt,vt){gt,n)) dH 1 

t 

(Q(E(v t )) - kt) ( < =^+ (guvt)) dH 1 . 
Here v t is the outer normal, tj the tangent, kt the curvature of r gt and it is the unique solution 
CE(v t ) : E(w)dz = - [ div n ({g t ,v t ) CE(v t )) -wdH 1 , Vtoe A{fl gt ) . 



/"(*) = is jr i9t>vt)=- I 2Q(E(v t ))dz+ I \d Tt (g u v t )\ 2 dU 1 



(4.1) 



r„ 



Remark 3.4 and Lemma 3.6 yield 

d 2 

Jfi^W' ut >\t=o ~ u ^V") "nyj ^ <-o\\\y, "i\\Hi{r h y 



/"(0) = —F{g u vt)\ t=0 = d 2 F(h,u)[g] > c ||<£, v)" 2 
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It will be shown later in Lemma 4.7 that, when S > is chosen to be small enough, the previous 
inequality implies 

(4-2) f\t) = ^H9uVt)>^\\{^t)\\l HVgt) forall te[0,l). 



dt 2 OT " " 2 

2tt 

•••• ''II'// 7 .,„,. Holds lor all / -I- II. I |. Mure |, 



It is now clear that \\{gt, Vt)\\ 2 H i(Y g j > c I Ifl 1 * 1 1 i 2 ([o 2tt)) holds for all i € [0,1]. Since J Q *g 2 d6 
J 27r /i 2 d# we have f *(h — g) 2 d0 = 2 J 27r h(h — g) d9 and therefore 



1 / r 2v i / r 2n \ 2N 

1 II- 2 jo 1 / / /, n2 



(4.3) || fe || W)) = ]jjOT -^ ^ (fc _ flr<W J J. 

Since f^ih — g) 2 d6 is very small we obtain from (4.3) that 

(4-4) ll5tlli=([o,2 OT )) > \\\K-9\\l*([ow 

From (4.2) and (4.4) we conclude that f"(t) > c\\h — g|| 2 2- Since (h,u) is a critical pair we 
have /'(0) = and therefore 

F(g,v)-F(h,u)=f(l)-f(0)= [ (l-t)f'(t)dt 



>c\\h- g\\ L 2 ([0 ^ }) / (l-t)dt 



- c \\h 112 

- 2 " #ll£ 2 ([0,27r))' 



which proves the claim when g is smooth 

J # 



When g £ Ch 1 '^) the claim follows by using a standard approximation. □ 



It remains to prove (4.2). The proof is based on a compactness argument and for that we have 
to study the continuity of the second variation formula (3.4). To control the boundary terms in 

(3.4) we need fractional Sobolev spaces whose definition and basic properties are recalled here. 
The function h is as in Proposition 4.1 and is its graph. 

Definition 4.2. For < s < 1 and 1 < p < oo we define the fractional Sobolev space W s,p (Th) 
as the set of those functions v € L p (Th) for which the Gagliardo seminorm is finite, i.e. 

(4.5) [»U = (J Th l h ' MHvWHz)) 1 " < oo. 
The fractional Sobolev norm is defined as \\v \\w=.p(r h ) '■= IMlLJ>(r h ) + M*,p;iV 

The space W~ s ' p (Th) is the dual space of W s ' p {Yh) and the dual norm of a function v is defined 

as 

IMIw- = .P(r h ) :=SUP {_^, ™ ffll ( Z ) I \\ U \\w°*r(dF h ) < 1 

We also use the notation H s (T h ) for W s ' 2 (T h ) for -1 < s < 1 and the convention W°' p (T h ) := 
L p (Th). By Jensen's inequality we have the following classical embedding theorem. 

Theorem 4.3. Let —1 < t < s < 1, q > p such that s — 1/p > t — 1/q. Then there is a constant 
C depending on t,s,p,q and on the C l -norm of h such that 

\\v\\ w *.«{T h ) < C\\v\\w»{T h y 
We also have the following trace theorem. 

Theorem 4.4. If p > 1 there exists a continuous linear operator T : W 1,p (Qh) — > VF 1_1 / p,p (rfc,) 
such that Tv = v\r h whenever v is continuous on The norm of T depends on C x -norra of h 
and 7. 

The next lemma will be used frequently. 
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Lemma 4.5. Let — 1 < s < 1 and suppose that v is a smooth function on I\ . Then the following 
hold. 

(i) If at C^Th) then 

\\av\\ W s, P(rh) < C||a||ci(r fc )||u||w".p(r fc ), 

where the constant C depends on p, s and the C x -norm of h. 

(ii) Ifty:Th^>- ^(Th) is a C 1 -diffeomorphism, then 

\\v° i Jt\\w:p(9(rh)) ^ c \\ v \\w^p(r h ), 
where the constant C depends on p, s and the C x -norms of h, \P and 

We will also need to control regularity of the elastic equilibrium. To this aim, the following 
elliptic estimate turns out to be useful, see [14, Lemma 4.1]. 

Lemma 4.6. Suppose (g,v) £ X(0) is such that 7 < g < R — 7, g e C^(R) and v € A(Q g ) 
satisfies 

(4.6) / CE(v) : E{w) dz = f : E{w) dz for every w € A(n g ), 

Jn g Jn g 

where f £ C (fi ff ; M 2x2 ). Then for any p > 2 we have the following estimate 

\\E(v)\\ wl „ (flt .flpx* ) + ||VCE(i;)|| j 
(4-7) ( v 

< C ( ||£J(«)|| Z ,2( f J g . M 2X2) + ||/||cl(n s ;M2x2) J , 

where T denotes the space of third order tensors and the constant C depends on 7, p and the 
C 2 -norm of g. 

We are now in position to give the proof of the inequality (4.2). To control the bulk energy 
we use techniques developed in [14]. The main difference is that we use directly elliptic regularity 
rather than deal with eigenvalues of compact operators. 

Lemma 4.7. Suppose that a critical pair (h,u) € X(uo) is a point of positive second variation 
with < h < Rq and = 1 • Then there exists S > such that for any admissible pair 

{g,v) € X Teg (uo) with II^Hl 2 = 1 an d \ \h — g\\c 2 (R) < 6 we have for 

h + t(g - h) 

9t ~ 



\h + t(g-h)\\ L 2 ' 
that 

(4-8) ^H9uV t )> C ^\\(g_t lVt )\\ 2 Hl(Tgt) for all te [0,1], 

where Vt is the elastic equilibrium associated to g t . The constant cq is from the Lemma 3. 6. 

Proof. Choose 7 > such that 7 < h < R — 7. Suppose that the claim is not true and there are 
pairs (g n ,v n ) S X rcg (u Q ) and t n e [0, 1] with 

\\h- g n \\c 2 (R) -> 

for which the claim doesn't hold. Denoting 

d ( h + t(g n -h) \ 8 2 ( h + t(g n -h) 

1 1 and g n := 1 



Jn - dt\ t=tn \\\h + t(g n -h)\\ L2 J * n ' dt 2 \ t ^\\\h + t(g n -h)\\ L2 
this implies 

(4-9) lim ll/ .- F " (g "' 2 t; " ) < ?. 

n ^°° \\\9ji,Vn)\\Hi( rg j 2 
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where by the notation T" {g n ,v n ) we mean 

F"{g n ,v n ) = - [ 2Q(E(v n ))dz+ I \d Tn {g n ,v n )\ 2 dH 1 



(d Vn Q(E(v n )) + k 2 n ) {gn,u n ) 2 dH 1 
( 4 ' 1Q ) + / (Q(E(v n )) - k n ) d Tn ((gn, v n )(gn, r n )) dH 1 

(Q(E(v n )) - k n ) [ ( .~ Vn) l + (g n ,v n )) dU 1 

\ \9n-> V n) — J 
= h+I 2 + h + h + h. 

Here v n is the outer normal to F Sn , r„ and k n are the tangent vector and the curvature of T g 
and v n is the unique solution to 

(4.11) / CE(v n ) : E(w) dz = - [ div T „ v n ) CE(v n )) ■ w dU 1 , Vw £ A(Sl gt ) . 



As in Proposition 3.1 we find C°°-diffeomorphisms : — > Q, gn such that : Th —> and 

ll*n ~ «rf|lc 2 (O h ;R2) < C\\h- g n \\ C 2 (R) . 

The goal is to examine the contribution of each term in (4.10) to the limit (4.9). We begin by 
proving that the contribution of I4 and I§ to (4.9) is zero. 

Notice that the C 2 -convergence of g n implies fc„ovI/ n — > k in L°°(rft,). Moreover, since v n solves 
the first two equations in (2.7) and sup n ||gn||c 2 ([o,27r)) < C, we have by a Schauder type estimate 
for Lame system, see [14], that there is a € (0, 1) such that 

(4.12) sup||w„|| c i, Q( a ;i . R 2 ) < 00, for Q' n = B iJo _ 7 \F 9n . 



\\E(uoy- 1 ) - ^( u „)||^ 1(P(ngn . M3> < 2) +||VC£( W o^- 1 )-VC J BK)|| ff _ 4(rgn;T) 



Next we prove the following elliptic estimate 
(4.13) 

<C\\h — 5n||c 2 (R)j 

where p > 2 and C depends on 7, p and the C 2 -norms of h and u. Indeed by the equations (2.7) 
satisfied by u and v n and a standard change of variables we obtain 

(4.14) f C(E(u o *" x ) - E(v n )) : E(w) dz = f f n : E(w) dz, G A(Sl gn ), 

f) 



where /„ € C 1 (fi s „;M 2x2 ) satisfies 



1 |c 1 (« £ ,„) < C||/l - 5n||c 2 (R) 

for C depending only on the C 2 -norm of u. Lemma 4.6 yields the estimate 

\\E(u o - £(« n )|| wl „ (neB . Ma „ ) + HVC^o - VCi?K)|| ff _i (r ^ ;T) 

< C (||£(uo l^ 1 ) - S(w„)|| L 2 (ngri . M 2x2 ) + \\h - g n \\ C 2 (R) ) . 
On the other hand, using w — u o ^J^ 1 — i>„ as a test function in (4.14), we obtain 

\\E(U O ) - £'(u„)|| I/ 2(Q gn . M 2x2) < C||/„|| i 2 (n9n . M 2x2). 

This concludes the proof of (4.13). 

By the trace theorem 4.4, Lemma 4.5 and (4.13) we obtain 

\\E(v n o *„) - E(u)\\ H i {rhM2x2) < C\\E(v n o $„) - B(u)|| Hl(nhiMa x a) 

< C||5n - /l llc 2 (R)- 

This estimate together with (4.12) implies v n o fy n — > u in C 1 '". In particular, we have that 

(Q{E(v n )) - k n ) o -> Q(B(«)) — k = X 
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uniformly, where A is a Lagrange multiplier. We may use this to estimate the term 7 5 in (4.10). 
By explicit calculations one easily obtains that ||(<7n) ^IIl 1 < C| I (Snj^n) 1 1 1^,2 and recalling that 
the functions g n and g n satisfy the volume constraint, as in Remark 3.4, 

{ r Un)2 x + (g n ,v n )dU l =Q 
we get 

(Q(E(v n )) - k n ) [ Vn \ + (g n ,u n )) dU 1 
\ {gn,v n ) — J 

(Q(E(v n )) - k n - A) f ( .~ Vn) l + (g n ,v n ) ) dU 1 

\ \9n,V n ) — J 

< C\\Q(E(v n )) -k n - ^L^(r gn )\\{9n,^n)\\ 2 L 2 (rgn y 
Using polar decomposition we have v n = ^5=s= and r„ = 9 ^J \ + g ,"f ~ ■ Since (g ni T n )(z) = 
7j^{9ni u n)(z) an d II jf^yllir 1 ^) < C we have as above that 

(Q(E(v n )) - k n ) d Tn ((gn, v n )(gn,T n )) dH 1 

< C\\Q(E(v n j) - k n - A|U« ( r 9 j||(4^„)||^ ( r g „)- 

Hence the contribution of the terms I4 and ^5 to the limit (4.9) is zero. 

The remaining terms 1\ , I2 and ^3 form a quadratic form. The goal is to show that 

(4.15) hm = hm ^ ^ 



n ^°° IK5n,^n>llffi(r 9 „) \\{±n^)WH^(T h ) 

where 

, 9ng-n 
ipn ~ 



h 

and Ujf, n solves 

(4.16) / CE(u^ n ):E(w)dz=- I div r ((V^, v) CE(u)) ■ wdH 1 , Ww e A(Q h ) . 

Notice that ipn satisfies the volume constraint J 27T hip n dO = then we may use the Lemma 3.6 to 
conclude that 

hm TTTi \TT2 - c o ' 

n-voo \\{lpn,ls)\\ Z Hl{rh) 

which then contradicts (4.9) and proves the claim. 

To show (4.15) we will compare the contribution of each terms in the quadratic form 

9 2 J(/i,u)y = - / 2Q(E{u iln ))dz+ f \d T {ip n ,v)\ 2 dH 1 - [ (dvQiEiu^ + k^i^rfdn 1 



with respect to the one given by Ii,l2 and ^3 in (4.10). 

We first point out that since g n — > h in C 2 we have that v n o <]/„ — >• v and r„ o vp„ — >• t in 
C 1 (r/ l ). Therefore from the definition of ip n we get 

u IKV^, v)\\m(r h ) _ 

™^°° IKSry v n)\\Hi(T gn ) 

and the convergence of li-, 

f \d T (^n,v)\ 2 dH 1 

lim -4^ = 1. 

\dr n (9n,Vn)? dU 1 
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The convergence of Ii follows from the equations (4.11) and (4.16). Indeed, by using standard 
change of variables, these equations yield 
(4.17) 

CE{u i>n of,- 1 ) -CE(v n j) : E{w)dz = [ (/„ E(u^ n o l^ 1 )) : E(w) dz + [ d n ■ w dH 1 

for any w G A(Q gn ). Here 

d n = div T «Vv», v)CE{u)) o f- 1 1^*,; 1 ! - div Tn ((gn, u n )CE(v n )) 
and /„ G L 2 (n gn ;M 2x2 ). For /„ we have 

(4-18) H/n||L~(n 3 „;M2x2 ) -> 0. 

By the estimate (4.13) we get 

Therefore, by Lemma 4.5, the choice of and from — l||c 1 (B) — > we have that 

Hr 9 „;K 2 )' 



9n 

(4-19) 11^1^-1(^^)11(^^)11^^) "►<>. 

Choose 

w(z) = (u^ o f ,- 1 - u„)(z) + Az + b 
as a test function in (4.17) where A is antisymmetric and b is a vector. This yields 
(4.20) 

Q{E(u i>n o f- 1 - v n ))dz 



+ 1 1 oL 1 1 i 1 1 w 1 1 i 

11 nH //-5(r B „;R 2 ) n M ff2(r g „;R 2 ) 



< C||/n||L~(n g , 1 :M 2 x 2 )||£ : (uVr 1 )lk 2 (nh;M 2 x 2 )||^(uV'r l ° *n _ ^n) 1 1 L 2 (0 Sll ;M 2x2 ) 
II "11^-2(1 

By the Theorem 4.4 we get that 



"'II ..i .. < C\\w\\ H X(Q gnm ^y 



'H2 (r 3ll ;K 2 ) 

As in the proof of Lemma 3.6 we choose A such that 

\\w\\Hi(n gn :M 2x2 ) ^ C \\ E ( w )\\L 2 (n gn ;M 2x2 ) = C\\E(u^ n o f" 1 - W„)||L 2 (f2 s „;M 2 x 2 ), 

by Korn's and Poincare's inequalities (choose b accordingly). The two previous inequalities and 
(4.20) yield 

WEiu^oi&n 1 -V n )\\ L 2(n gn . M 2x2) <C'^||/„|| L oo( naji . M 2 X 2)||£:(M V , n )|| L 2( nh . M 2x2) + | | d n \ \ R _ 1 _ R2 J . 

Arguing as in (3.20) we may estimate 

ll^( u V>„)IU 2 (n h ;M 2 x 2 ) < C \\{-tj)n,v)\\ H -L^ h y 

Therefore using (4.18) and (4.19) we deduce that 

II^OvJIIja^M 2 * 2 ) ~ I |-EQ>n) 1 1^2(^^2x2) 
I Kffn^n) 1 1 ffi(r 8n ) 

This proves the convergence of 1%. 

We are left with the term J3 in (4.10). We need to show that 

{d Vn Q{E{v n ))+k 2 n ) {g^v n ) 2 dU x - f {d v Q{E(u))+k 2 ) {^v) 2 dH 

Due to the C 2 -convergence of g n and the C 1 -convergence of ^ we just need to show 
(4-21) \\d Un Q(E(v n )) o * n _ d v Q{E(u))\\ H _ h[Th) 0. 



\\(9n^n)\\ H i(r g J -> 
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This will be done as [14, Proposition 4.5]. For every tp £ H^iT^) we have 

= S h { CE (£r) °^- CE (£) ) : wo *•> * 

+ ^ CE [-^) : (E(v n )o^ n -E(u))^d'H 1 

< \\(vcE( Vn )) o * n - vc^(«)|| ff _j (rh;T) ||(^K) o * n ) <p\\ H h {rhM2x2)) 

+ C\\E(v n ) o - S(u)|| i 2 (r(i . M 2x2 ) ||( / 3|| L2(rh) 

where the constant depends on C 2 -norms of u and h. Fix p > 2. By the definition of Gagliardo 
seminoma, Holder's inequality, Theorem 4.3 and Theorem 4.4, we obtain 

\\(E(v„)oH! n )(p\\ i 

< C\\(E(v n ) o * n )\\ L «, {rh . U 2 X2)) \\<p\\ H i (rh) + o *n)|l w ^ ,^ 2x lblli»(r h ) 



< C M| o *„)|| LTC(rh . M 2x2 )) + \\(E(v n ) o 



2p 



w p- 2 (n h ;M 2 >< 2 )/ "H 2 (r h ) 
By repeating the previous argument for we obtain by (4.12) and (4.13) that 
\\VQ{E(v n )) o - VQ(S(«))|| H _i (rfciEa) 

( 4 - 22 ) < C (H(VC^K)) o tf n _ VC£?(u)|| ff _ 4 (rh . T) + \\E{v n ) o * n - S( w )|| L 2 (rh;M 2x2 ) ) 
< C||fli„ - /i||c 2 (R)- 

Since v n o fy n —> v in C 1 , (4.22) implies (4.21). This concludes the convergence of the term I3 and 
completes the proof. □ 

5. Local minimality 

This section is devoted to prove the main result of the paper, the local minimality criterion. 
Namely, we show that if a critical point (h, u) € X ies (uo) has positive second variation, then it is 
a strict local minimizer in the Hausdorff distance of sets and a quantitative estimate in terms of 
the measure of the symmetric difference between the minimum and a competitor holds. Due to 
the sharp quantitative isoperimetric inequality, the exponent 2 in (5.1) is optimal. 

Theorem 5.1. Suppose that (h,u) € A reg (ito) is a critical pair for T with < h < Rq. If the 
second variation of T is positive at Qi,u), then there is 6 > such that for any (g,v) € X(uq) 
with |O s | = \Q h \ and < du^g u ^g^h) < S it holds 

(5.1) T(g,v)>T(h,u) + c\n g An h \ 2 , 

for some c > 0. 

The proof is based on a contradiction argument and follows some ideas contained in [14], [9] 
and [1]. Assume that (h n ,u n ) is a contradicting sequence, i.e., 

HhmO <F{h,u) + c \n hn An h \ 2 and < d n (T hn U Z hn ,T h ) < -. 

n 

The idea is to replace (h n ,u n ) by (g n ,v n ), which minimize an auxiliary constrained-penalized 
problem, and to prove that (g n ,v n ) are sufficiently regular to apply the C 1, ^minimality criterion 
to get a contradiction. As auxiliary problem we choose 

min I F(g, v) + A| \Q g \ - \Q h \\ + yj (|n fl An fc | - s n ) 2 + e n : {g,v) e X(u ) , g<h+^\, 

where the second penalization term will provide the quantitative estimate in (5.1) and the obstacle 
g < h + l/n plays a key role in proving the regularity of (g n , v n ). 



18 



GIUSEPPE MARIA CAPRIANI, VESA JULIN, AND GIOVANNI PISANTE 



The regularity proof is divided in three steps. In Lemma 5.6 we prove that g n is Lipschitz using 
some geometrical arguments. Then, in Lemma 5.7, we show that g n is a quasiminimizer for the 
area functional which in turns implies C 1 '"-regularity. Finally, we deduce the C ' -regularity in 
Lemma 5.8, by using the Euler-Lagrange equation for (g n ,v n ). 

The following isoperimetric-type result will be used frequently in this section. The proof can 
be found in [1, Lemma 4.1]. 

Lemma 5.2. 

(i) Let f £ C|?(]R) be non-negative and let g £ BV#(JL), then there exists a constant C, 
depending only on f , such that 

(ii) Suppose D is a set of finite perimeter. Then 

P(D U BJx)) - P(BJx)) > -\D\ , 

r 

where P stands for the perimeter. 

We will also need the following property of concave functions. 

Lemma 5.3. Suppose that f n £ C 1 (IR) and f £ C 1 (M.) are such that ||/„ — /||l<»(h) 0. If f n 
are concave then it holds 

Wfn-f\\ci c (M) ->0. 
Proof. Let R > and fix e > . Since f £ C 1 (K) we find S > such that 

f(S + x)-f(x)>f(x)6~e6 
for every |je| < R. On the other hand, since f n are concave, we have 

fn{S + x) - f n (x) 
g < /nW- 

Hence 

nx) - m < MzM+^MiAW) + £ < 2e , 

when n is large. The reverse inequality f' n {x) — fix) < 2s follows from a similar argument. □ 

We begin the study of the properties of solutions of the auxiliary problem by proving an exterior 
ball condition. 

Theorem 5.4. Let h £ CT(R) such that < h < R and r < min{l/(A + 1), I/Hk/JI^}, where 
Kh is the curvature ofTh- Then for every c, e £ [0, 1] and n £ N any solution of the problem 

(5.2) min|j-( 5! u)+A||O s | - \Q h \\ + c^(\Sl g AQ h \ - e) 2 + e : (g,v) £ X(u ), g<h+^ , 

satisfies the uniform exterior ball condition, i.e., there exists r > such that for every z £ dF g 
we find zq for which B t {zq) C R 2 \ F g and dB r (zo) n (T g U S fl ) = {z}. 



Proof. The proof follows the argument from [14, Lemma 6.7]. Recall that dF g = T g U S s . Given 
a ball B r (zo) define the half circle S r (zo) — dB r (zo) n {z £ M 2 : (z — zq,zo) < 0}. Assume that 
there exists a ball B r (zo) C K 2 \ F g such that S r (zo) intersects T g U S g in two different points 
z\ = (pi,0i) and z 2 = [p2^i)- When r < l/W^hW^ it is clear that the arc / = f(0) of S r (z ) 
connecting z\ and z^ satisfies f(9) < h(0) for 8 £ (01,82). Let g be defined by g = f for 8 £ (8\, 82) 
and g — g otherwise. Denote by / the arc of T g U E s connecting z\ and Z2 and by D the region 
enclosed by / U /, see Figure 5.1. 
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Figure 5.1. 




Notice that 



|fi ff Afi h | - e) 2 + e - \/(\n g An h \ - ef + e 

(\n s An h \- £ ) 2 -{\n g An h \-£) 



(5. 3) ^/(\n g An h \-e)2+e + ^/(\n g An h \ - ef + e 

\n g An h \ + \n g An h \ - 2e){\n § An h \ - \n g An h \) 



< 

\\n g An h \ -e\ + \ \n g An k 

< l%An fl |. 



Since tt g A£l g — D and fl g C fi ff we see that 



T(g,v) + A\\n g \-\n h \\+cJ(\Q g AQ h \ 

(5.4) 



< F(g, v) + A| |n 9 | + |n h || + <y (|n 9 Afi h | - e) 2 + e + HHf) - n.f) + (A + 1)|D| ■ 

Since P(D U B r (z )) - P(B r (z )) < U l (f) - H(f), Lemma 5.2 implies 

HHf)-H(f)<~\D\. 

Hence, since r < 1/(A+1), the inequality (5.4) contradicts the minimality of (g, v). The conclusion 
now follows arguing as [8, Lemma 2] or [12, Proposition 3.3, Step 2]. □ 

Lemma 5.5. Let h,c,e and n be as in the previous theorem. Suppose (g,v) € X(uq) is any 
minimizer of (5.2). Then there exists A$ > ; independent of c,e and n, such that if A > Ao then 

\n B \>\a h \. 

Proof. We argue by contradiction supposing that |fi fl | < |fi/i|. We observe that there exists 
< r < 1 such that, if we define fl g = Br \ rF g , we have = \Slh\- Moreover, since 

r 2 r 27T 



we get 

= f nR 2 -\n h \ 

1 \7rR 2 -\n g \ 



< i. 
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Clearly £l r g — Q 3r for g r (9) = rg(9). Define the function v r : Vt gr — > R 2 as 

'wo(^o) if rRo < \z\ < Ro 
if) if g r (^) <\z\<rR . 



v r (z) 



Since Q, gr D il g , we see that |f2 9r .Af2 ff | = |f2^| — \flg\- Using the inequality (5.3) we have, for A 
sufficiently large, that 



F{g r ,v r )+A \\n gr \ - \n h \\+cJ(\n gr ACl h \-E) 2 +e 



- T{g, v) - A||Q 9 | - p h \\ - c^(\n g AQ h \ - e) 2 + e 

< [ Q{E(v r ))dz-A(\n h \-\n g \) + c\n gr An g \ 

JrR <\z\<R 

<C(l-r)-(A-l)(|fi h |-|n,l) 

<c(\n h \-\n g \)-(A-i)(\n h \-\n g \)<o, 

which contradicts the minimality of {g, v). □ 

In the following we study convergence properties of solutions for constrained obstacle problem 
(5.2). 

Lemma 5.6. Let h be as in Theorem 5.4- Assume g n G BV#(R) is such that g n < h + 1/n and 
it satisfies the uniform exterior ball condition. If 

(5.5) g n -> h mL 1 and lim •H 1 (r„ USJ= H 1 {T h ) , 

n— >oo 

then g n — > h in L°° . Moreover, for n sufficiently large, g n are uniformly Lipschitz continuous. 
Proof. Here we follow an argument from [14, Theorem 6.9, Steps 1 and 2]. We claim that 

sup \g„ — h\ — > as n — > +oo . 

Let us first note that T gn U S ffn is a connected compact set. Up to a subsequence, we can assume 
that T gn U converges in the Hausdorff distance to some compact connected set K. The 
continuity of h and condition (5.5) imply that T/j C K. By Golab's semicontinuity Theorem (see, 
e.g. [3, Theorem 4.4.17]) and assumption (5.5) we see that 

T-^i^h) < T-O^iK) < lim n 1 (T gn U^ gn )=n 1 (T h ). 

n— >+oo 

This implies that 'H 1 (if \T^) = 0. Since K is connected, it follows from a density lower bound (see, 
e.g. [3, Lemma 4.4.5]) that K = Th- Now the claim follows from the definition of the Hausdorff 
metric and from the continuity of h. 

From the previous claim we see that, for n sufficiently large, it holds 7 < g n < Rq — 7 for 
some 7 > small, independent from n. Hence, since the polar coordinates mapping is a C°°-local 
diffeomorphism far from the origin, the representation in polar coordinates of F 9n (still denoted 
by F 3n ) satisfies the uniform exterior ball condition up to changing the radius r to f € (0, 1) by a 
factor depending only on 7. Now we prove that g n are L-Lipschitz with L < S ||/i||c 1 (R)- 

We argue by contradiction and assume that there exists 9 and 9k — > 9 such that 

, llm a\ - ~ C 1 ^) 

fc^oo \9k — 0] r v ' 

and set z — (9,g n (9)). Without loss of generality we may assume that 9k > 9. By the uniform 
exterior ball condition we find a ball B?(zo) C R 2 \ F 3n such that dBf(zo) n (F gri US 9n ) = {z} and 

( M 1 \ 4 

Z0 = z+ "{vTTm>vTTm)> for M ^fiwic MR ) 
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Let z' e dBf(z ) such that 

' VM 2 - 3 2 



zo - r 



yi + A/ 2 ' Vl + M 2 ^ 
We write z' =: z + f (u>i, W2) with 



M - VM 2 - 3 , -1 

w\ = ; — > and W2 = = < 

and since Bf(zo) cl 2 \ F 9n we have g n (9 + fwi) < g n (9) + rw^. Setting S n = sup R |/i — g n \ and 
recalling ||/i||c l (R) < M/4 we get 

fM fM 
/i(0 + ru>i) > - —v>i > g n (8) - S n - — Wl . 



Therefore we deduce 



' M 

h(9 + fwi) - g n (9 + fwi) > -S n - r ( -j-toi - w 2 

f ( M 



= - i " + 7TTm( 1 -T( M -^ 



-3 



r / 3M \ 

= —S n H . 1 . > S n 

VI + M 2 \ 4(M + VAf 2 - 3) J 
which, for sufficiently large n, is a contradiction. □ 

In the next lemma we show C '"-regularity of the minimizer for the penalized obstacle problem. 

Lemma 5.7. Let h be as in Theorem 5.4 and (g n ,v n ) € X(uo) be any minimizer of the problem 
(5-6) 

mhx(j 7 (g,v)+A\\n g \ - \Q h \ \ + Cy J(\n g An h \ - e n ) 2 + e n : (g,v) e X(u ) , g < h+ , 
where c € [0, 1] and e n — > 0. Assume also that g n — > h in L 1 and that 

lim U^iYa US, ) =H 1 (T h ) and sup/ Q(E(v n ))dz < +00. 

Then for all a £ (0, |) and /or n Zarc/e enough g n £ C 1,a (M.), the sequence {Vv„} is equibounded 
in C a > a (n gn ;M 2x2 ), and g n h in C^ a {R). 

Proof. From Lemma 5.6 we infer that g n is sufficiently regular to ensure a decay estimate for Vi>„. 
Indeed, for z Q G T gn there exist c„ > 0, a radius r n and an exponent a n £ (0, 1/2) such that 



I, 



\Vv n \ 2 <c n r 1+2a " 



'B r (z )nn 9 „ 

for every r < r n . This follows from the fact that v n minimizes the elastic energy in Q gn and the 
boundary T gn is Lipschitz, see Theorem 3.13 in [12]. 

Since g n is Lipschitz, we may extend v n in B r (zo) such that 

(5.7) / |VS„| 2 <c„r 1+2a ", 

JB r (z ) 

where v n stands for the extension. 

For r < r n , denote by z' r and z" the two points on T gn n dB r (z ) such that the open sub-arcs 
of T gn with end points z' r , z and z", zq are contained in r 9n n dB r (z Q ). Setting z' r = g n (6' r )a(6' r ) 
and z" = g n {9")o-(9"), denote by I the line segment joining z' r and z" and define 



9n(<?) 



ffnW fl€[0,27r)\(^,^) 
mm{h(6) + 1,1(6)} 9 e 

where 1(9) is the polar representation of I. 



22 



GIUSEPPE MARIA CAPRIANI, VESA JULIN, AND GIOVANNI PISANTE 



By (5.7) and by the minimality of the pair (g n ,v n ) we have 
(5.8) W x (r 9 „ n B r (z )) - V. 1 (Tg n n B r (z )) < C n r 1+2a < 

Indeed we can estimate 

o>^{g n ,v n ) - Hgn,v n ) + A(\\n 9n \ - \n h \\ ~ \\n Sn \ - \n h \\) 



+ c \J(\Q gn An h \ - s n ) 2 + e n - yJ(\Sl Sn Ml h \ - e n ) 2 + e 

>U l {T gn r\B r {z Q ))-'H\T gn r\B r {z Q ))- f Q{E{i n ))dz-{K + l)TTT 2 

JB r (z ) 

>U l {T 9n n B r (z )) - U l {T- gn n B r (z Q )) - C n r 1+2a - 
We will show later that 

(5.9) 7~L 1 {Tg n n B r (z )) — H X (l) < Cr 2 . 
Now the inequality (5.9) together with (5.8) gives us 

(5.10) u 1 (r gn n B r (z )) - n\i) < Cr 1+2a - 

and the desired C 1,Q! -regularity follows from a classical result for quasiminimizers of area functional 
(see Theorem 1 in [20]) once we observe that 

% X (0 — ml B r (z )) : FA0 9 „ g B r (z )} . 

The proof of (5.9) is a consequence of the C 2 -regularity of h and goes as follows (see Figure 5.2): 



< 1 [ r (~g 2 n -l 2 )de+- [ r (~g' n + l')(~g' n -l')d9 
7 h' T 7 Je' r 

< -\B r (z )\ + C [ r (~g' n -l')d6, 
7 7 J9' r 

where C depends on the Lipschitz norm of g n and I in the interval (0' r ,6") and 7 is a positive 
constant with 7 < ming h. 
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To estimate the last term we first note that either the set {h + l/n < 1} is empty or there exists 
0o G (9^,9") such that g' n (8o) — 1'(9q) = and using second order Taylor expansion around 9q we 
easily get 

\g' n -l')d9<Cr 2 

where C depends on the C 2 -norm of h. 

Now we claim that g n converges to h in C 1 -norm. As in the proof of Lemma 5.6 we will work in 
the plane (6,p) and we recall that the subgraph of g n , still denoted by F 9ri , satisfies the uniform 
exterior ball condition. From the C^-regularity and the uniform Lipschitz estimate, in the Lemma 
5.6, we obtain sup n H^nHc 1 < °°- Hence, from the uniform exterior ball condition we conclude 
that at every point there exists a parabola touching g n from above. In other words, there is C > 
such that for every 9 it holds for P{9) = g n (9 ) + g , n (9 )(9 - 9 ) + C (0 - 9 ) 2 that 

min(F- 5 „) = (P-g n )(9 Q ) = 0. 



This implies that g n are uniformly semiconcave, i.e., for every n the function 

9^g n (9)-C9 2 

is convex. We may now use the Lemma 5.3 to conclude the desired C 1 -convergence of g n . 

The convergence of g n to h in C 1 -norm allows us to use a blow-up method (see [14, Theo- 
rem 6.10]) to infer the uniform estimate 

(5.11) / \Vv n \ 2 < c r 1+2 ° 

JB r (z ) 

for any a G (1/2, 1) and for all r < r$ where cq and r$ are independent of n. 

Once we have (5.11), we can repeat the argument used to prove (5.10), replacing (5.7) by (5.11), 
to infer 

H\T gn n B r (z )) - H\l) < Cr 1+2a . 
This implies a uniform estimate for the C 1,a -norms of g n for a G (0,1/2) (see for instance [9, 
Proposition 2.2]). The (^'"-convergence of g n now follows by a compactness argument. 

To conclude the proof we have just to observe that, since v n is a solution of the Lame system 
in flg n , we can apply the elliptic estimates provided in [14, Proposition 8.9] to deduce that Vv n is 
uniformly bounded in C°' a (Ti g7i ,R 2 x R 2 ) for all a E (0, 1/2). □ 

Lemma 5.8. Let (h, u) G X reg (uo) be a critical point of T such that < h < Rq, and let (g n , v n ) 
be as in the previous lemma. Suppose that Vw„ — Vu weakly in L^ oc (r2/ l ; IR 2 x M 2 ), 

lim / Q(E(v n )dz = / Q(E{u))dz and lim ^ g " A ^ = 1. 

n ~^°°Jn gn Jn h n ^°° £n 

Then g n G C 1,1 (IR) and g n — > h in C » (M), for n sufficiently large. 

Proof. From Lemma 5.7 we know that g n — > h in C 1 ' a (R). Therefore for large n there exist 
diffeomorphisms $„ : Cl gn — > fj/, such that — > id in C 1,a . Let Br be any ball of radius 
R G (Rq — maxR h, Rq). Since 

sup {|K|| C l,a ( n_ )} < 00 

by the convergence Vv n — 1 Vm we have that 

(5.12) Vv n o^- 1 ^Vu in C°' a (n h nB R ; M 2x2 ). 

To prove the claim denote /„ := {9 G [0, 2ir] \ g n (9) < h(9) + 1/n =: h n (9)}. Since I n is open, 
we may write J„ = {j^Li(a^, &"). Notice that 

(5.13) g' n (9) = h' n {9) = h'{9) on [0, 2tt] \ I n . 

If /„ is empty, the claim is trivial. Therefore we may assume that /„ ^ 0. Since g n G C 1,Q (E), we 
can write the Eulcr-Lagrange equation for (g n ,v n ) in the weak sense: 

(5.14) k gn (6) = Q(E(v n ))(6,g n (6))+(3 n (6,g n (e)) + \ n , 9 G /„. 
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Here 



y/{\n gn An h \-e n ) 2 + e n 



sign (xn h - Xn„ n ) 



and A„ is some Lagrange multiplier. Notice that from the assumption limn-^ £ n 1 |Sl flrl Af7^ | = 1 
it follows that 



(5.15) 



< A- 



in„ An h 



o. 



y/(\n gri ASl h \-e n ) 2 +e r . 
Recall the Euler-Lagrange equation for (h, u) 
(5.16) k h (6) = Q(E(u)){8,h(9)) + X OD . 

We will show that X n — > Aoo . Notice that for the curvature in polar coordinates it holds 

91 + 2 9n - 9n9n 



■ fin 



fin 



fin 



(gl + g'?)^ "" \Vgl + 9nJ Val + fin 

Hence, multiplying (5.14) by g n , integrating over I n and using (5.16) yield 



Q(E(v n )(e,g n (0)) + p n (9,g n (6)) + \ n \g n d6 = ^ k gn g n d6 



fin 



de 



E 

oc 

E 



Vgn+g'nJ Vdn + 9n 

9n(b?) g'JaV 



V9 2 n m + g' n 2 m V9 2 n (a?) + g' n 2 (K)J 
KM) _ h' n (a?) 



fin 



a™ V 9n + g'n 



de 



fin 



kh h n d9 



de 



- I d'i 

'«? V 9n + 9n 



fin 



Q{E{u){0,h{0j)+*°° hd0+ I (k hn h n -k h h)+ _ 

Ir, V 9n + 9n V + 



dB. 



Recall that h n = h + l/n. Therefore by (5.12), (5.15) and the previous calculations we conclude 
that 

lim — - / \ n g n - Xoohde = , 

n-+oo \l n \ J In 

which clearly implies A„ — > A^. 

From (5.13) and (5.14) we conclude that g n £ C 1,1 ^). Moreover by the equations (5.13), (5.14) 
and (5.16) together with the convergences (5.12), (5.15) and A„ — > we conclude that 



kn — > kh 



This in turn gives us the convergence 



in L°°. 



in L° 



a 



Now we are in position to prove the main theorem of this section. 

Proof of Theorem 5.1. Step 1: We show first that (h, u) is a strict local minimizer, i.e., we prove 
the claim without the estimate on the right-hand side of (5.1). 

Observe that from the results of the previous section we may assume that (h, u) is a C ,1 - 
local minimizer. The result will follow once we prove that C^-local minimality implies the local 
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minimality. Arguing by contradiction let us assume that for any n € N there exist (h n , u n ) G X(uq) 
with \Clh n \ = \^h\ such that 

T{h n .u n ) < F(h,u) and < d n {T hn U E 9 „, T h ) < -. 

Consider the sequence {g n , v n ) € X(uq) of minimizers of the following penalized obstacle problem 

1 



mmU : (g,v)+k\\Q. a \-\Q, h \\ : (g,v) e X(u ), g < h+ ^ 

for some large A. Since (h ni u n ) and (h,u) are clearly competitors, we have that 

F{gn,V n ) < F(h n , u n ) < F{h, u) . 

By the contradiction assumption we may assume that (h ni u n ) ^ {h,u). 

By the compactness property of X(uq) there exists (g,v) such that, up to subsequences, 
(gmVn) (g,v) in X(uo). Let (f,w) € X(uq) with / < h, by the lower semicontinuity of T 
and the minimality of (<?„,«„), we get 

H9,v) + A||fi fl | - \il h \\ < liminf [j"(5„,u„) + A||J2 9 „| - \il h 

(5.17) L 

<r(f,w) + A\\si f \-\n h \\. 

Choosing (/, w) = (h, v) in the previous inequality, we obtain that 

(5.18) H 1 (T g ) + A\\n g \-\n h \\ <H\T h ) 

When A is sufficiently large, (5.18) and Lemma 5.2 imply that g — h. Moreover, we observe that 
from (5.17) it follows that (h,v) minimizes J- in the class of all (f,w) € X(uq) with f < h. In 
particular v must coincide with the elastic equilibrium u. 
Choosing (/, w) — (h,u) in (5.17), we deduce 

lim H 1 (r 9 „UE s J=W 1 (r /l ), 

n — ^oo 

lim I Q{E(v n ))dz= ( Q(E{u))dz. 

From Lemma 5.8 we get g n — > h in C 1,:L (R). 

We only need to modify g n such that it satisfies the volume constraint. We simply define 
9n{0) '■= gn{9) + 5n where 5 n are chosen so that = |f2/»|. By Lemma 5.5 it holds \£l Bn \ > \£lh\ 

and therefore 6 n > and Q gn C ^ 9 „- Hence v n is well defined in £l gri and (g n , v n ) is an admissible 
pair. 

Since h > and g n — > h uniformly, we have g n > 7 for some 7 > independent from n and 
5 n — > 0. We may estimate 



Jo 

1 f 27r 

< - / 2 5t A + ^ dfl 
7 Jo 



and 



Therefore whenever A > - we have 
— 7 



27T 1 /.27T 

(ft, + 8 n f -9 2 n de=-J^ 2g n S n + 5 2 n d6. 



(5-19) U l {T~ gn ) -U l {T gn ) < A||fi fi J - |fi 9 J|. 

The claim now follows, since by the choice of g n and by (5.19) we have 
^(g n ,v n ) = F(gn,v n ) + A||%J - |f2 h || 

<-F(0n,«n)+A||n ft J-|fifc|| <T(h n ,u n ) <F{h,u). 
This contradicts the fact that (h, u) is a strict C^-local minimizer. 
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Step 2: We will now prove the theorem. The proof is very similar to the first step. Arguing by 
contradiction we assume that there are (h n ,u n ) € X(uo) with |f2/t„| = \£lh\ such that 

Hhn,u n ) <F{h lU ) + c \n hn An h \ 2 and < d n (T hn U £ g „, T h ) < - . 

Denote e n := \Q hn A£l h \. Notice that dfi(T hn U T, gn ,T h ) —> implies Xn hn —> Xn h m L 1 and 
therefore e n — ^ 0. 

This time we replace the contradicting sequence (h n ,u n ) by {g n ,v n ) G X(uq) which minimizes 

mm(j : (g,v) + A\\n g \ - \Sl h \\ + yJ(\n g AQ h \ - e n f + e n : (g, v) £ X(u ), g < h + i 

By compactness we may assume that, up to a subsequence, (g n ,V n ) —¥ (g,v) in X(uq). By 
completely similar argument as in Step 1 we conclude that (g, v) = (h, u) whenever A is sufficiently 
large. Moreover, we have that 

lim H 1 (T gn UE gn )=n 1 (T h ), 

n—too 

lim / Q(E(v n ))dz= I Q(E{u))dz. 
To conclude that g n — > h in C 1,1 (IR), it remains to show that 
(5.20) lim I"'-*"*' = 1 

n— >oo £ n 

and apply Lemma 5.8. 

Suppose that (5.20) were false. Then there exists c > such that |f2 ffn AO^| — e n \ > ce n . Use 
the minimality of (g n ,v n ) and the contradiction assumption for (h n ,u n ), we obtain 



F{g n ,v n )+K\\VL gn \- \n h \\ + (\n gn An h \ - e n ) 2 + e n 

( 5 - 21 ) < «n) + 

< J"(/i, u) + c e^ + 

Now we observe that from [13, Proposition 6.1], for A sufficiently large, (h, u) is also a minimizer 
of the penalized problem 

F(g,v)+A\\n g \-\Sl h \\. 

Hence we have 

(5.22) F(h,u) < F(g n ,v n ) + A\\Q gn \ - \Q h \\. 

Combining (5.21) and (5.22) we get 



V 7 c 2 e 2 n + e n < y (\Q gn AQ h \ - e n ) 2 + s n < c a e 2 n + y/s^, 

which is a contradiction since e n — > proving (5.20). 
Arguing as in (5.21) and by using (5.20) we obtain 

F{g n ,v n ) + A\fl g J\ - p h \ \ < F{h n ,u n ) + Vi^- ^(\n 9n An h \ - e n ) 2 + e n 
(5-23) < F{h 

, u) + c e n 

<T{h,u) + 2c \n gn An h \ 2 , 

when n is large. 

As in Step 1 define g n (0) :— g n (6) + S n where 6 n > are such that \£l gn \ — \Clh\- By choosing 
A large enough we have 

A 
2 

Therefore since 



(5.24) ^(TsJ-n 1 ^) < -\\n- gn \-\n gn \ 



n gn AQ h \ 2 < 2\n Sn An h \ 2 + 2\n gn An g j 2 = 2\n gn An h \ 2 + 2\\n g j - \n h 
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we have by (5.23) and (5.24) that 

F(Mn,v n ) < F{g n ,v n ) + — \\Slg n \ - 

< F{h,u) + 2c lO^Afi^l 2 - -||fi ff J - \Vt h 



< F(h,u)+4eo\Sl §n ASl h \ 2 - |||J)„„| • • <>,,| ,~4<„|i<>., I - \<h 

<T(h,u) + Ac \ng n An h \ 2 , 



,2 



when n is sufhciently large. This contradicts Proposition 4.1 when c$ is chosen to be small 
enough. □ 

6. The case of the disk 

In this section we consider the particular case when a radial stretching is applied to a material 
with round cavity F = B r . We prove that the disk remains stable under small radial stretching. 
This result is similar to the case of flat configuration in [14]. The main difference to the flat 
case, where the minimal shape is a rectangle, is that the curvature of the disk is nonzero and 
therefore the second variation formula becomes considerably more complicated. Instead of trying 
to explicitly write the second variation, we use fine estimates to find a range of stability. 

The Dirichlet boundary condition has the form of radial stretching, 

(6.1) u Q {pa{6)) = aRoa(d) for p > R , 

where a 6 R is some constant. The region occupied by the elastic material is the annulus 
A(Ro,r) := Bji \ B r . For uq as above we say that (h, u) 6 X(uq) is a round configuration if 
h(9) = r and u is the elastic equilibrium associated to h. 
For the next theorem we define 



m : = i + i 



p + A t 2 

ITri 

Recall also the definition of the ellipticity constant 77 = min{/i, p + A}. 
Theorem 6.1. Let 

r :=sup(<<i? | (l + ^logY^ ) > V 



t ) ~ 4/i 



and define the function G : M — > [— 00, Rq) as 



G(a) := sup it < Ro | Hog ( ^ ) p 2 {t) > '' 



t ) r K "' ~ 32(p + X) 2 a 2 
If r G (ro, Ro) and a € K satisfy 
(6.2) r > G{a) , 

then the round configuration is a strict local minimizer of J- under the volume constraint. 

The elastic equilibrium u can be explicitly calculated. Indeed, because of the symmetry we can 
write 

u(pa(9)) = f(p)a(9) 
and applying the first equation in (2.7) we have 

r(p) + />)_M = o. 
p p z 



This can be easily solved 



f{p) = -+bp, 
P 
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for some a.kl. To find a and b observe that 
(6.3) 

(f'(p)cosH+lf S m 2 9 (f(p)- £f)sin9cose\ ( ^Wj 

C ( " } ^(/'(^-^sinflcos* f(p)sin^ + mcos^J + H /(P) + ^Ao 1 

Therefore, the second equation in (2.7) gives 

(2 M + A)/'(r) + A^ = 0. 

r 

This and the Dirichlct condition (6.1) yield 
(6.4) 4 = ^6 and b= a 



P P{r) ' 

It is trivial to check that the round configuration is a critical point of J- '. To prove the Theorem 
6.1 we need to show that the round configuration is a point of positive second variation. To this 
aim, let us explicitly write the quadratic form (3.16). By (6.3) and (6.4) we have 

„„, , .,, , , / sin 2 9 — sin#cos#\ .,, , > 
CE(u) = ibtfi + A) . a a 2 „ = 4b(fi + X) t<E)t, 

w y ' \— sin 9 cos 9 cos z 9 J v ' 

on the boundary dB r . Hence, for ip € iTiQR.), we have 

div r ( i/) C£(u) ) = 46(^ + A) (- i/) v + d T (±, i/) r) 

and the equation (3.17) for it^, becomes 

(6.5) / CE(u^):E{w)dz = -4b(n + X) I (-(?/>, i/)(ui, i/) + (d T {ip, v))(w, r)) tfH 1 . 

JA(R,r) JdB r 

Moreover, in the case of a round configuration the elastic energy is 

a 2 

(6.6) Q(E(u)) = 2(p + X)b 2 + 2n — 

P 

and therefore, by (6.4), we get 



du Q(E(u)) = - 8iti + X)2 - on dB r 



Hence, (3.16) becomes 



( 6 - 7 ) , M2.2 



d 2 F(h,uM=- 2Q(E(u^,))dz + \d T (±,v)\ 2 cM 1 

JA{R,r) JdB,. 

( 8(p + \) 2 b 2 _ 1_ 
as,. V M r r 2 



where € A(Br \ B r ) solves (6.5), and ^ satisfies ij}d9 = 0. 

Now the goal of is to prove that (6.7) is positive whenever the assumptions of Theorem 6.1 are 
satisfied. The main obstacle is to bound the first term in (6.7) which will be done by using the 
equation (6.5). To this aim we need the following simple lemma, which we prove to keep track of 
the optimal constant. 

Lemma 6.2. Suppose that v £ W 1 ' 2 (A(R 0l r); M 2 ) is a continuous map with v = on <9-Br and 
A is a matrix. Then for w(z) = v{z) + Az we have that 



[ IwfdH 1 <rlog(^ / 
JdB r \ r J J a 



A(Ro,r) 



Dv A 

Ro-r 



dz. 
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Proof. Consider w in polar coordinates. Fix an angle 9 and integrate over [r, R ] 



ARq<t{6) — w(ra(9)) = / Dw (pa (9)) a (6) dp 



Ro 



which implies 



Ro 



\w(ra(9))\ < / Dv(a(9)) 



Ro-r 



-A 



dp. 



Integrate over 9 and use Holder's inequality to obtain 



\w(pa(9))\ 2 d9 < 



2tt / r R 



\ Jr 
2tt / r R< 



Dv(pa(9)) 



-A 



Rn-r 



dp d9 



< 



lo. 



1 

o \Jt P 

Rq 



dp ■ 



Dv(pa(9)) - A 

Ro-r 



pdp) d9 



A(Ro,r) 
2tt 



Dv 



Rn-r 



-A 



dz. 



The inequality follows from f gB \w\ 2 dW 1 = r f Q n \w(r,9)\ 2 d9 



□ 



Proof of Theorem 6.1. As we stated before, by the local minimality criterion it is enough to prove 
that the second variation of T at (h, u) is positive. Suppose that ip e i/i(R) satisfies %p d9 = 
and ip 0. Without loss of generality we may assume ip to be smooth. To estimate the first term 
in (6.7) we claim that 



(6.8) 2 



Q(E(u^))dz< 32(/i + A)2b2 rlog 

A(R Q ,r) V 



Rq 



(i>,vf + \d T (i>,v)\ 2 dU X 



To this aim, choose w(z) ~ u^(z) + Az as a test function in (6.5) where A is antisymmetric, to 
obtain 



(6.9) 



A(Ro,r) 



Q(E(u^,))dz = -Ab(n + \) / (-{jfr, v)(w, v) + d T {±, v)(w, r)) dU 1 



dB r 



<4b(p + X)[ I (±,v) 2 + \d T (±,v)\ 2 d% 1 

JdB T 



1/2 



IwfdH 1 



8B r 



1/2 



Apply Lemma 6.2 to w to get 



(6.10) 



< rlog 



dB r 



Rq 
r 



A(Ro,r) 



Duj, - — A 

R — r 



dz. 



Let R k — > oo and for every k choose an antisymmetric A k such that 



>A(R k ,r) 

By Theorem 2.5 we get 

r 

Du^p - — A k 

lA(R k ,r) H—r 

Together with (6.10) this yields 



Du^ - — A k dz = 

R — r 



A(R k ,r) 



Dui + A k dz. 

v R — r 



dz<C k \E( Ui ,)\ 2 dz = C k \E(uip)\ 2 dz. 



A{R k ,r) 



A(R,,r) 



[ \wfdn 1 <rlog(^-)c k / \E(u^)\ 2 dz. 

JdB r \ r J JA(R ,r) 



Since C k — >■ 4 as R k — >■ oo we have that 

IwfdH 1 < — lo 



2^1^ ^^{Ro 



c)B r 



V 



A{Ro,r) 



Q(E(u^))dz. 



Now (6.8) follows from (6.9) and from the previous inequality. 
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We estimate (6.7) by using (6.8) and obtain 

d 2 F(h, uM > -Wfj- 1 ^ + \fb 2 r log (^j J {£, v) 2 + \d T {±, v) | 2 dU 1 

\d T (M\ 2 dH 1 + f ( 8 ^ + x ) 2 l_ 1 \ it)U f dH i 

JdB r \ A 1 r r J 



(6.11) = / \d T {±,u)\ 2 -\{±,v) 2 dn x 

JdB r r 

-32rr 1 (M + A) 2 6 2 rlog (^j J \d T {±,v)\ 2 dU 1 

Let us first treat the last term in (6.11). For every r > ro we have that 

d 2 F{h,u)[i>]> (^l-32? ? - 1 ( M + A) 2 fe 2 rlog^^^ \d T {±,v)\ 2 - ^{±,v) 2 dU 1 . 
Furthermore, if (6.2) is satisfied, then 

1 - 32?T 1 0u + A) 2 6 2 r log \ >0. 

By the definition (3.1) we see that (ip, v) — ip ( <7 ~ 1 (|f|))" Hence, by the Wirtinger's inequality, 
we get 

/ \d T {±,v)\ 2 - ~{±,v) 2 dU 1 = 1 W{6)\ 2 -\m\ 2 d9>Q. 

JdB r r r JO 

which concludes the proof. □ 

At the end of the section we study the global minimality of the round configuration. We begin 
with the following remark. 

Remark 6.3. Suppose that Rq and ro are as in Theorem 6.1 and fix a g R and a small e > 0. 
Then for every re [ro + e, Ro] such that r > G(a) + e the proof above actually gives 

d 2 T(h,uM >c 2 W{9)\ 2 -\^{9)\ 2 de + Cl \ij(9)\ 2 d9, 



for some small < c\ < C2, independent of r. Using the Wirtinger's inequality we get 

d 2 F(h,uM >C ||V>||lri ([ o,2,r)): 

for Co depending only on R , r , a and e. This is a uniform version of the Lemma 3.6. 

We can use this uniform bound of the constant Co to prove a uniform local C^-minimality of 
the round configuration for r G [ro+£, Ro] with r > G(a)+e. Indeed, arguing as in the Proposition 
4.1 and in the Lemma 4.7 we conclude that there is 5 > such that for any (g, v) E X(uq) with 
\F g \ = \B r \ and \\g - r|| c i,i( R ) < S it holds 

where u r stands for the elastic equilibrium associated to the disk B r . 

The previous remark enables us to prove the global minimality of the disk when the volume of 
the annulus is small. 

Proposition 6.4. Suppose that Ro is the radius of the large ball and uo is the Dirichlet boundary 
conditions as in (6.1) with fixed a > 0. There exists r g i b < Ro such that for every r € (r g i b, Ro) 
the round configuration, with a disk B r , is a global minimizer of J- under the volume constraint. 
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Proof. We argue by contradiction and assume that there exist a sequence of radii r n /*■ Rq and 
a sequence (k n ,W n ) £ X(uq) of minimizers of F under the volume constraint |ftjt n | = |^4(i?0j ?"n)| 
such that 

F(k n ,w n ) < F(r 

where u n stands for the elastic equilibrium relative to v n . Since (k ni w n ) minimizes F we im- 
mediately have that H 1 (T^ ri U — » 2hRq. Therefore, since Fk„ is connected, we deduce that 
£n '■= d-H(Tk n U Sfc n ,r rii ) — > as n — > oo. 
We may calculate the elastic equilibrium 

u n (p, 9) = + b n f^j cr(0), 

where 

, fi,( 1 + Xr l\ A V + A 2 r 

o n = H « and a„ = r n b n . 

By the Remark 6.3 we have that for large n it holds 

d 2 F(r n ,u n )[i/;} > co||^||Hi(aB r „)! 

for Jg 2 ^ ip d8 = 0, where cq is independent of n. 

We note that u n is also the elastic equilibrium in the annulus A(R,r n ), for any R > Rq, with 
respect to its own boundary conditions on dB R , v(R,8) = u n (R,6). For R > R we define 

F R (g,v) = I Q(E(v))dz + H 1 (T g )+2H 1 (X g ) 

JB R \F g 

and 

X R {u n ) = {(g,v) | g e W # (M), « € fl^,, (M 2 \ F g ;M 2 ), u = u n outside B^} . 
Consider the estimate (6.11) for d 2 F Rl (r n ,u n ) [tp], i.e., replace Ro,r and b by Ri,r n and 6 ra . By 
continuity we may choose R\ close to i?o such that 

d 2 F Rl (r n ,u n )[^] > yllV'll^cas^), 

for f Q 2,r ijjdO = 0. Arguing as in the Remark 6.3 we conclude that (r n , u n ) is a local C 1,1 -minimizer 
of F Rl uniformly in n, i.e., there is 6 > 0, independent of n, such that for any (g,v) € X Rl (u n ), 
with ||gf — r„|| c i,i(R) < S, it holds 

(6.12) F Rl (g,v) > F Rl (r n ,u n ). 

Define 

w n (z) if zeB Ro \F kn 
u n (z) if z € i?o) ■ 



w„(z) := 



By the assumption on (fc„, io„) it holds 

(6.13) F Rl (k n ,w n ) < Fri {r n ,u„). 
Suppose that {g n ,v n ) is a solution of the problem 

min{ F Rl (g, v) + A||F fl | - \B r J\ : (g,v) S X Rl (u n ), g < r n +£„}, 

where A is large. Arguing as in Lemma 5.6, Lemma 5.7 and Lemma 5.8 we conclude that g n — » i?o 
in C M (IR). In particular, \\g n - r n || c i.i (R ) -> 0. 

By the minimality of (g n , u n ) we have that F Rl (g n , v n )+A\ \F g \ — \B Tn | | < F Rl (k n ,w n ). Defining 
fjn = gn + S n such that |i 7 'g n | = \B Tn \ we obtain, as in (5.19), that 

(6.14) F Rl (g n , v„) < F Rl (g n ,v n ) + A\\F g \ - \B Tn \ \ < F Rl (k n , w n ), 

when A is large enough. Moreover S n — > 0. Hence \\g n — r n \\cui(m ~^ an d therefore (6.12), 
(6.13) and (6.14) imply 

F Rl (r„, u n ) < F Rl (g n , v n ) < F Rl (fc„, w n ) < F Rl (r„, u n ), 

which is a contradiction. 
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